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I. INTRODUCTION 



The Q-state Potts model p] in two and three dimensions exhibits a rich variety of critical 
behavior and is very fertile ground for the analytical and numerical investigation of first- 
and second-order phase transitions. With the exception of the two-dimensional Q = 2 Potts 
(Ising) model in the absence of an external magnetic field, exact solutions for arbitrary Q are 
not known. However, some exact results at the critical temperature have been established 
for the two-dimensional Q-state Potts model for no magnetic field. From the duality relation 
the ferromagnetic critical temperature is known to be T c = J / k B \\i{l + ^fQ) for the isotropic 
square lattice. The ferromagnetic Potts model has a first-order phase transition for Q > 4 
in two dimensions ||, and has been the most important means to understand first-order 
transitions which are omnipresent in nature, not only in condensed-matter physics but also 
in subatomic physics, for example, the deconfining or chiral symmetry breaking transition in 
quantum chromodynamics[§, |]] and phase transitions in the very early Universe ||. In three 
dimensions the ferromagnetic Q-state Potts model for Q > Q c k, 2.45(1)0 has a first-order 
transition, and the three-state Potts model is related to the order parameter (the expectation 
value of the Wilson line) for the deconfinement transition in quantum chromodynamics0. 

On the other hand, the antiferromagnetic Potts model is much less well understood than 
the ferromagnetic model. One of the most interesting properties of the antiferromagnetic 
model is that for Q > 2 the ground-state is highly degenerate and the ground-state entropy is 
nonzero. The Baxter conjecture^] for the critical temperature of the Potts antiferromagnet 



(J < 0) on the square lattice, T c = J/fceh^v^ — Q — 1), gives the known exact value for 
Q — 2, a critical point at zero temperature for Q = 3, and no critical point for Q > 3. In three 
dimensions the Potts model on the simple-cubic lattice has the antiferromagnetic ordered 
phase for Q < 40. It is now commonly accepted that the order- disorder phase transition of 
the three-state Potts antiferromagnet on the simple-cubic lattice belongs to the universality 
class of the three-dimensional XY model || [K|. Recently Rosengren et al. [Jll, [1^] claimed the 
existence of two additional phase transitions below the order-disorder phase transition for the 



simple-cubic three-state Potts antiferromagnet. However, Kolesik and Suzuki|fl3|, Heilmann 
et a/.|14j], and Oshikawa||15|| reported different results in disagreement with observations by 
Rosengren's group. 

By introducing the concept of the zeros of the partition function in the complex magnetic- 
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field plane, Yang and LeeflllJ proposed a mechanism for the occurrence of phase transitions 
in the thermodynamic limit and yielded a new insight into the unsolved problem of the Ising 



model in an arbitrary nonzero external magnetic field. It has been shown [TO, |T7| , [I8| , |T9| , 
2"0| , |2TJ that the distribution of the zeros of a model determines its critical behavior. Lee 



and Yang|T7| also formulated the celebrated circle theorem which states that the partition 
function zeros of the Ising ferromagnet lie on the unit circle in the complex magnetic-field 
(X = e^ H ) plane. In three dimensions the properties of the zeros in the complex X plane 
of the Ising model have been studied by exact enumeration [22] and series expansion |p3 



However, for the Q-state Potts model with Q > 2 the partition function zeros in the complex 
X plane lie close to, but not on, the unit circle with the two exceptions of the critical point 



X = 1 (H = 0) itself and the zeros in the limit T = 0f2|, ||, g§. 

Fisher p7j emphasized that the partition function zeros in the complex temperature 



(Y 



plane are also very useful in understanding phase transitions, and showed that 



for the square lattice Ising model in the absence of an external magnetic field the zeros in 
the complex Y plane lie on two circles (the ferromagnetic circle Yfm = — 1 + V2e td and 
the antiferromagnetic circle Yafm = 1 + V2e ld ) in the thermodynamic limit. In particular, 
using the Fisher zeros both the ferromagnetic phase and the antiferromagnetic phase can 
be considered at the same time. Fisher also showed that the logarithmically infinite spe- 
cific heat singularity of the Ising model results from the properties of the density of zeros. 
The zeros in the complex Y plane of the simple-cubic Ising model have been studied exten- 
sively to understand the analytical properties and critical behavior of the model by exact 



enumeration |28|, |29| , |30| , pl| |3^[ and to determine the critical point and exponents by Monte 
Carlo simulation [pi] , [33| , |34|, [35], [36|]. Recently it has been shown that for self-dual boundary 
conditions near the ferromagnetic critical point Y c = 1/(1 + y/Q) the zeros of the Q-state 
Potts model on a finite square lattice in the absence of an external magnetic field lie on the 
circle with center —1/(Q— 1) and radius y/Q/{Q — 1) in the complex Y plane, while the an- 
tiferromagnetic circle of the Ising model completely disappears for Q > 2 



38, M 
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Using the distribution of the zeros in the complex temperature plane of the square-lattice 
Potts model, the critical behavior in both the ferromagnetic and antiferromagnetic phases 
has been studied |39|, 51 1, and the unknown properties of the model for nonzero external 
magnetic field have been invetigated|]39|]. The partition function zeros in the complex Y 
plane of the three-state Potts model on the simple-cubic lattice have been studied to under- 
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stand the ferromagnetic critical behavior by exact enumeration for the 3x3x9 simple-cubic 
lattice and by Monte Carlo simulation pl|, |43[ . 



In this paper we study the partition function zeros of the Q-state Potts model for 2 < Q < 
7 on the simple-cubic lattice to unveil some of the rich structures of the model. The partition 
function zeros of the three-dimensional Potts model have never been obtained for Q > 3 in 
the literature. In the next section we describe briefly the microcanonical transfer matrix to 
evaluate the density of states, from which the exact partition function of the simple-cubic 
Potts model is obtained. In Section 3 we discuss the ferromagnetic and antiferromagnetic 
phase transitions and the complex-temperature singularity (nonphysical critical point) of the 
simple-cubic Potts model using the partition function zeros in the absence of an external 
magnetic field. We also discuss the properties of the complex-temperature singularity using 
the low-temperature series expansion. In section 4 we study the unknown properties of the 
simple-cubic Potts model for nonzero magnetic field using the partition function zeros and 
the high-field, low-temperature series expansion. We discuss the Yang-Lee edge singularity 
and the Fisher edge singularity, which are nonphysical critical points for nonzero magnetic 
field, and relationship between the Yang-Lee and Fisher edge singularities. 

II. DENSITY OF STATES 

The Q-state Potts model in an external magnetic field H q on a lattice G with N s sites 
and Nf, bonds is defined by the Hamiltonian 

N s 

7~{-Q — — J ^2 5((Ti, (Tj) — H q ^2 Si^ki <Z)j (1) 
(id) k=1 

where J is the coupling constant (ferromagnetic model for J > and antiferromagnetic 
model for J < 0), indicates a sum over nearest-neighbor pairs, 5 is the Kronecker 

delta, <Ji = 1,2, ...,Q, and q is a fixed integer between 1 and Q. The partition function of 
the model is 

Z Q = E e"^, (2) 

{<T„} 

where {<J n } denotes a sum over Q Ns possible spin configurations and (3 = (ksT)" 1 . If 
we define the restricted density of states with energy < E < Nb and magnetization 
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< M < N s by 

Oq(E, M)=J2$(E-J2i 1 - <rj))WM ~ £ ?))> ( 3 ) 
which takes on only integer values, then the partition function can be written as 

Z Q (Y,X) = Y- N » £ £ Q Q (E,M)Y E X M , (4) 

E=0 M=0 

where Y = e _/3J and X = e f3Hq , and states with E — (E — Nf,) correspond to the 
ferromagnetic (antiferromagnetic) ground states. From Eq. (4) it is clear that Zq(Y,X) is 
simply a polynomial in Y and X. If we sum the restricted density of states over M, we 
obtain the density of states 

n Q (E)= Y.V Q (E,M), (5) 



M=0 

whose sum over E is equal to the Q Ns spin configurations, 

N b 

E ^q(E) = Q Ns - (6) 

E=0 

In the absence of an external magnetic field, the partition function of the model is given by 

N b 

Z Q {Y) = Y~ Nb £ Q Q (E)Y E , (7) 

E=0 

which is again a polynomial in Y. 



Here we describe briefly the microcanonical transfer matrix (/xTM) p2| , ^T], 0, |45], ^7 
to obtain exact integer values for the density of states Qq(E) of the Q-state Potts model 
on an L x L x iV simple-cubic lattice with periodic boundary conditions in the x-y plane 
(horizontal area of L x L) and free boundaries in the z direction (height of N). First, an 
array, u^, which is indexed by the energy E and variables <jj, 1 < i < L 2 for the first x-y 
plane of sites is initialized as 

LxL 



(1) 



[E- ax, (7 2 , a LxL ) = 8(E - E [1 - S(at, a i+1 )]). (8) 

8=1 

Now each spin in the plane is traced over in turn, introducing a new spin variable from the 
next x-y plane, 

^ 2 \E-a[,a 2 ,...,a LxL )=Y / ^ (1 HE-[l-6(a[,a 1 )] ] a 1 ,a 2 ,...,a LxL ). (9) 



This procedure is repeated until all the spins in the first plane have been traced over, leaving 
a new function of the LxL spins in the second plane. The bonds connecting the spins in 
the second x-y plane are then taken into account by shifting the energy, 

LxL 



J 2 \E; a[, a' 2 , a' LxL ) = Q&>(E - £ [1 - 5(4, a' i+1 )}; a[, a> 2 , a' LxL ). (10) 

i=i 



This procedure is then applied to each x-y plane in turn until the final (iVth) plane is 
reached. The density of states is then given by 

n«(£) = EE- E ^ N \E;a[,a 2 ,...,a' LxL ). (11) 

The permutation symmetry of the Q-state Potts model allows us to freeze the last spin 
°lxl — 1 of each x-y plane. Now we need to consider only Q LxL ~ 1 possible spin configura- 
tions in each plane instead of Q LxL configurations, and we save a great amount of memory 
and CPU time. For example, Table I shows the density of states Qq=-?(E) of the seven-state 
Potts model on the 3x3x3 simple-cubic lattice. 

The above procedure can be modified in a straightforward way to calculate exact integer 
values for the restricted density of states Q(E, M)p5j [47]]. For the first x-y plane one has, 
instead of Eq. (8), 



LxL LxL 



u {1) (E, M- a x , a 2 , o~lxl) = 5(E-J2l 1 - o i+1 )])5{M - £ 6{a k , q)). (12) 

i=l k=l 

The spins in the first plane are traced over in sequence just as in Eq. (9) leading to a new 
function of the spins in the second x-y plane; 

Cu (2 \E, M- a[,a 2 , a LxL ) = £ u {1) {E - [1 - 6{a[, at)}, M - 5(a[, q); a x , a 2 , a LxL ). (13) 

Once all the spins in an x-y plane have been traced over, the function for the next plane is 
completed as in Eq. (10) 

u^(E, M; a[, a' 2 , a' LxL ) = G>&>(E - E [1 - 5{a[, a' l+1 )}, M; a[, a' 2 , a' LxL ), (14) 

i=i 

and the restricted density of states is given by 

n Q (E,M) = J2J2- E ^ N \E,M;a[,a' 2 ,...,a' LxL ). (15) 
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III. PARTITION FUNCTION ZEROS IN THE COMPLEX TEMPERATURE 
PLANE 



In this section we discuss the properties of the Q-state Potts models on the simple- 
cubic lattice using the partition function zeros {^o} or {^0} of the models in the complex 
temperature (Y = e _/3J or A = e@ J ) plane which are obtained by solving the equation 

Z Q (Y ) = Y - N > £ n Q (E)Y E = (16) 

E=0 

for the exact partition functions on finite lattices. Figure 1 shows the partition function 
zeros in the complex Y = e _/3J plane of the Q-state Potts models for Q=2, 4, and 7 on 
the 3x3x3 simple-cubic lattice with periodic boundary conditions in the x-y plane and 
free boundary conditions in the z direction. In Figure 1, near the positive real axis, the 
partition function zeros show one-dimensional distributions for all the Potts models which 
become denser as Q increases. Since the value of latent heat of a system is proportional 



to the density of zeros in the complex temperature plane |19[, for Q > 3 the latent heat 
increases at the transition temperature as Q increases, and the first-order phase transition 
becomes stronger. Moreover, the first zeros of these linear distributions will cut the real 
axis and be the ferromagnetic (FM) critical points in the thermodynamic limit. In Figure 
1, as Q increases, the first zeros become closer to the real axis and to the values of the FM 
critical temperatures estimated in the literature (Table II). This fact suggests that the first 
zero obtained for a smaller size of simple-cubic lattice approach very closely the real critical 
point in the limit Q — > 00 where the mean-field theory is known to be exact [[K], [KJ . 

Even for the three-dimensional Ising model the critical temperature has never been known 
exactly except for some conjectures, for example, the Rosengren conjecture pDJ. However, 



for the simple-cubic lattice Hajdukovic conjectured|y]] the FM critical temperature 



of the Q-state Potts model based on two assumptions. The first assumption is that Yf M (Q) 
reduces to the mean-field result Yf M = Q^ 1 ^ 3 in the limit Q — > 00 148|, f|9|]. And the second 
assumption is that in d dimensions Yf M (Q) is the solution of an equation 

Y~ d - Cl (d)Y- d ' 2 - c 2 (d) =Q-2, (18) 
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where c\{d) and 02(d) have to be determined, for example, c\ = 2 and C2 = 1 for d = 2. As 
shown in Table II the Hajdukovic conjecture is close to numerical estimates in the literature 
for 2 < Q < 7. However, for Q = 1 the Hajdukovic conjecture, Y FM = 0.7937..., completely 
disagrees with the numerical estimate 0.753[52", [53]. We can guess a simpler form than the 
Hajdukovic conjecture for the critical temperature, 

Y ' FM(Q) = WTW- (19) 

satisfying the first assumption. Eq. (19) reasonably agrees with numerical estimates for 
1 < Q < 7 as shown in Table II. Even though the Hajdukovic conjecture and Eq. (19) 
are incorrect expressions for the FM critical temperature of the Q-state Potts model on the 
simple-cubic lattice, we can infer that the true expression for the critical temperature may 
be something between the Hajdukovic conjecture and Eq. (19). 

Recently using cluster- variation calculations Rosengren and Lapinskas[l l[] claimed the 
existence of two additional phase transitions at Trli = 0.787(— J/ks) and Trl2 = 
0.7715(— J/ks) below the well-known antiferromagnetic (AFM) critical temperature 
t afm = 1 2259(7) (-J/k B ) p| for the three-state AFM Potts model (J < 0) on the simple- 
cubic lattice. Using the Monte Carlo method Kundrotas, Rosengren, and LapinskasJT^ 



obtained T^rl — 0.68(— J/ks) and supported their earlier claim. However, Kolesik 
and Suzuki(|Tj|] and Heilmann et a/. [13] observed different results by Monte Carlo simu- 
lations which are in disagreement with observations by Rosengren's group, and obtained 
Tkh — l-0(— J/ks) as the intermediate transition temperature in the region below T^ FM . 



Furthermore, based on a renormalization-group picture Oshikawa||15|| argued that there ex- 
ists only one transition at T^ FM for the three-state AFM Potts model on the simple-cubic 
lattice. Figure 2 shows the partition function zeros in the complex A = e^ 3 plane of the 
three-state Potts model on the 3 x 3 x 30 simple-cubic lattice with periodic boundary con- 
ditions in the x-y plane and free boundary conditions in the z direction. The interval 
0<A<1(0<T< 00) is antiferromagnetic (J < 0) while the interval 1 < A < 00 
(00 > T > 0) ferromagnetic (J > 0). There exists a clear one-dimensional distribution of 
zeros near Akrl = 0.2298, in Figure 2, which favors the existence of additional phase tran- 
sitions below A^ FM = 0.4424. However, we do not know if the first zero of the distribution 
near Akrl cuts the real axis in the thermodynamic limit, and if there exist the true phase 
transitions below A^ FM . 
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We have also obtained two pairs of the complex-temperature (CT) singularities (Fig- 
ure 2) using the low-temperature series 0, |62], |63| for the three-state Potts model on the 
simple-cubic lattice. Figure 2 exhibits accumulations of partition function zeros near the CT 
singularities. These accumulations imply that some thermodynamic quantities may diverge 
at the CT singularities which are not physical. The critical exponents associated with the 
CT singularities, which we call the CT critical exponents, are defined in the usual way, 

C° CT ~(l-—) , (20) 



l CT 

"CT 



Yct ■ 



}~r) ' <21) 
and 

/ Y \ ~~1CT 

Xor~(l-vH ' < 22 > 



Yct ■ 

where Yct is the location of a CT singularity, and C CT , m CT and \ct are the singular parts 
of the specific heat, spontaneous magnetization, and susceptibility, respectively. Here, the 
superscript represents no external magnetic field. Because the density of zeros near a CT 
singularity is given by[f27j 

9(Y) ~ (l - y-) , (23) 
the density of zeros diverges at the CT singularity if acT > 1- 



Guttmann et al. |64|, |65|, £6] studied the locations of the CT singularities and the CT 
critical exponents for the Ising model in detail, introducing a new method of low-temperature 
series analysis. For the Ising model on the simple cubic, body-centered cubic, and face- 
centered cubic lattices there exist the CT singularities, closer to the origin than the physical 
critical points, which control the low-temperature series and spoil the analysis of the real 
critical behavior. Guttmann et al. |M| reported that for the Ising model on the simple-cubic 
lattice the CT singularity is located at (Yct) 2 = -0.2860(2) (\(Y CT ) 2 \ < (Y C FM ) 2 = 0.4120) 
with cict ~ 1? ~\ < Pct < ~JE' an d Ict ~ §■ Itzykson et a/.|p0fl studied the zeros 
of the exact partition function in the complex temperature plane of the Ising model on 
the 4x4x4 simple-cubic lattice, and found that the distribution of zeros supports the 
real existence of the CT singularities. By the analysis of the longer series they obtained 
acT = 1.25(15), Pct = —0.05, and 7ct = 1.1(1) for the Ising model on the simple-cubic 
lattice, and argued that the values of the CT critical exponents of the Ising model are 
same for the simple cubic, body-centered cubic, and face-centered cubic lattices and satisfy 



acr + 2(3ct + 7ct — 2. Finally, Itzykson et al. conjectured that the behavior at the CT 
singularities is universal. Recently Guttmann and EntingpTfl obtained (Yct) 2 = —0.2853(3) 



with acr ~ 1.03, Pct ~ —0.01, and ^ct ~ 101 using the low-temperature series generated 
by the finite-lattice method of series expansion for the simple-cubic Ising model. They 
proposed occt = 1, Per = 0, and ^ct = 1 with logarithmic corrections, satisfying the 
Rushbrooke scaling law acr + 2Pct + Jct = 2. 



Using Dlog Pade approximants|>8[ to the low-temperature series []58[ of the specific heat, 
spontaneous magnetization, and susceptibility for the three-state Potts model on the simple- 
cubic lattice, we found the locations of two pairs of the CT singularities 

= 0.0254(11) ± 0.4948(9)i (24) 

and 

= -0.4981(45) ± 0.2453(35)z. (25) 

Since \Y$\ = 0.4955(9) and \Y$\ = 0.5552(42), these CT singularities lie closer to the 
origin than the estimated values of the FM critical point shown in Table II. The CT critical 
exponents for the simple-cubic three-state Potts model are same at Yqj, and Yqj, as shown 
in Table III, and satisfy the scaling law acr + ^Pct + Ict ~ 2. Because in three dimensions 
the values of the CT critical exponents for the three-state Potts model are very close to those 
for the Ising model, it is difficult not to suggest that the values of the CT critical exponents 
are independent of Q in three dimensions. That is, our results support the universality of 



the CT singularities in three dimensions! 5C]. For the two-dimensional Q-state Potts model 



on the square lattice it has been shown [59] that two pairs of the CT singularities of the 



model are mapped each other through the dual transformation 



which determines the FM critical temperature 

= rrw (27) 

Therefore, one expects the true expression for the FM critical points of the simple-cubic 
Potts model to provide the mapping between Y^ and Yqj,. 
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IV. PARTITION FUNCTION ZEROS FOR AN EXTERNAL MAGNETIC FIELD 



In this section we discuss the properties of the simple-cubic Q-state Potts models for 
nonzero magnetic field using the high-field, low-temperature series expansions [7C, 71, 72, 
[73], |74| and the zeros of the exact partition function Zq(Y,X) on a finite lattice in the 
complex temperature (Y = e _/3J ) or magnetic-field (X = e^ Hq ) plane. Figure 3 shows the 
partition function zeros in the complex X plane of the three-state Potts model on the 3x3x3 
simple-cubic lattice with periodic boundary conditions in the x-y plane and free boundary 
conditions in the z direction. The zeros lie close to the unit circle at lower temperatures, 
but they move away from the positive real axis and the unit circle as the temperature is 
increased. The zeros at T = are uniformly distributed on the circle with radius (Q — \y/ Ns 
which approaches unity in the thermodynamic limit, independent of Q, whereas the zeros 
at T = oo are iV s -degenerate at X — 1 — Q, independent of lattice size[|4], |5|, |26| . 



The edge zero, which we call the first zero, and its complex conjugate of a circular 
distribution of zeros in the complex X plane cut the positive real axis at the physical critical 
point X c = 1 for T < T c in the thermodynamic limit. However, for T > T c the first zero 
does not cut the positive real axis in the thermodynamic limit, that is, there is a gap in the 
distribution of zeros around the positive real axis. Within this gap, the free energy is analytic 



and there is no phase transition. Kortman and Griffiths [g^] carried out the first systematic 
investigation of the magnetization at the first zero, based on the high-field, high-temperature 
series expansion for the Ising model on the square lattice and the diamond lattice. They 
found that above T c the magnetization at the first zero diverges for the square lattice and is 
singular for the diamond lattice. For T > T c we rename the first zero as the Yang-Lee edge 
singularity. The divergence of the magnetization at the Yang-Lee edge singularity means 
the divergence of the density of zeros, which does not occur at a physical critical point. 
Fisher [ff5| proposed the idea that the Yang-Lee edge singularity can be thought of as a 
new second-order phase transition with associated critical exponents and the Yang-Lee edge 
singularity can be considered as a conventional critical point. The critical point of the Yang- 
Lee edge singularity is associated with a 3 theory, different from the usual critical point 
associated with the 4 theory. The crossover dimension of the Yang-Lee edge singularity is 
d c = 6. The edge critical exponent o e = l/5 e = (d — 2 + r] e )/(d + 2 — rj e ) at a Yang-Lee edge 
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singularity X e for Y > Y c is defined by 

m ; ~ (i - 1)-, ,28) 

where m e is the singular part of the magnetization in the complex X plane. The study 
of the Yang-Lee edge singularity has been extended to the classical n- vector model|7B|, the 
quantum Heisenberg modelJTB], the spherical model[]77|], the quantum one-dimensional trans- 



verse Ising model JFg], the hierarchical model |79|, the one-dimensional Potts model [gC], |8l 



branched polymers|]82[|, fluid models with repulsive-core interactions |83|, |84j], etc. Dhar|85 



calculated the edge critical exponent a e in two dimensions by solving a particular model of 
three-dimensional directed animals and mapping the solution to the hard hexagon model. 



Using Fisher's idea and conformal field theory, Cardy |36j studied the Yang-Lee edge sin- 
gularity for a two-dimensional 3 theory. It is generally accepted that the value of the edge 
critical exponent a e depends only on dimension. The values of a e are known to be a e = — \ 
in one dimension and a e = — ||85|, |86| in two dimensions, whereas a e ~ 0.08[]83| in three 
dimensions. 

The partition function zeros in the complex X plane for real Y have been extensively 
studied and well understood. However, the partition function zeros in the complex Y plane 
for real X are much less well understood than the zeros in the complex X plane. Figure 4 
shows the partition function zeros in the complex Y plane of the three-state Potts model for 
an external magnetic field X > 1 (H q > 0) on the 3x3x3 simple-cubic lattice with periodic 
boundary conditions in the x-y plane and free boundary conditions in the z direction. For 
X = 1 some zeros distribute along the one-dimensional locus near Y c , and the first zero of 
the distribution cuts the positive real axis at Y c in the thermodynamic limit as explained 
in the previous section. As H q increases, all the zeros move away from the origin. In the 
limit H q — > oo (X — > oo) the field H q favors the state q for every site and the Q-state Potts 
model is transformed into the one-state model whose the partition function zeros are the 
roots of the equation Y~ Nb = 0[[3£j. That is, \Y\ for all the zeros increases without bound 
as X increases. Note that for X > 1 there is accumulation of the zeros as we approach the 
first zero. Matveev and Shrock|87j studied the properties of the thermodynamic functions 
at the first zero in the complex temperature plane for the square-lattice Ising model with 
X > 1. They observed that the thermodynamic functions and the density of zeros diverge 
at the first zero when X > 1, and that the values of the critical exponents associated with 



12 



the first zero for X > 1 are nearly independent of X and satisfy the Rushbrooke scaling 
law approximately. We call the first zero in the complex temperature plane for X > 1 as 
the Fisher edge singularity The edge critical exponents a e , (3 e , and 7 e at the Fisher edge 
singularities are given by 

C e ~(l--) , (29) 



Y 



m Y e ~ (1 



YJ 



\/3e 

(30) 



Y \ _ > 

Xe-il-TT) , (31) 



and 

v. ~ I 1 — 

Y e 

where Y e is the location of a Fisher edge singularity, and C e , rrv^ and \ e are the singular 
parts of the specific heat, magnetization, and susceptibility, respectively, for X > 1. Kim and 
Creswickp9|1 extended the study of the Fisher edge singularity to the Q-state Potts model 
(Q > 3) on the square lattice. They found that the values of the edge critical exponents at 
the Fisher edge singularities are nearly independent of Q (Q > 2). 

In three dimensions to study the critical behavior at the Fisher edge singularity we have 
used the high-field, low-temperature series expansions for the Ising model 170, [n], |72| and the 



three-state Potts model [73, 73 1 on the simple-cubic lattice. Table IV shows the locations (Y e ) 
and the edge critical exponents (a e , f3 e , and j e ) of the Fisher edge singularities estimated 
by Dlog Pade approximants to the specific heat, magnetization, and susceptibility series for 
x = 100, 200, and 500. The edge critical exponents, which are independent of X, have the 
values of a e = 1.0 ~ 1.1, f3 e ~ —0.4, and 7 e ~ 1.0 for Q = 2 and for Q = 3 similar values of 
a e and j e and f3 e ~ —0.6. So in three dimensions we seem to have a strong violation of the 
scaling law a e + 2(3 e + 7 e = 2. Another interesting observation is that there exist apparently 
three pairs of the Fisher edge singularities in three dimensions (Figure 4 and Table IV). The 
estimates of the edge critical exponents at each of the three edge singularities indicate that 
they are equal within error bars. 

On the other hand, everything can be evaluated exactly in one dimension, and there exist 
only one pair of the Fisher edge singularities characterized by the edge critical exponents 
a e = |, (3 e = — |, and 7 e = |, all independent of Q. These values satisfy the scaling law 
a e + 2/3 e + 7 e = 2. Table V shows the locations and the edge critical exponents of the Fisher 
edge singularities for the Q-state Potts models on the square lattice. In two dimensions we 
have also found two pairs of the Fisher edge singularities which have the same values of the 
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edge critical exponents within error bars. The edge critical exponents have the values of 
a e = 1.13 ~ 1.24, P e = —0.11 ~ —0.20, and 7 e = 1.11 ~ 1.26. In two dimensions the edge 
critical exponents also satisfy the scaling law a e + 2(3 e + 7 e = 2. Now we can conclude that 
the number of pairs of the Fisher edge singularities is equal to the dimensionality. 

Until now the Yang-Lee and Fisher edge singularities have been studied and understood 
independently. However, the Yang-Lee edge singularity seems to have a direct relationship 
to the Fisher edge singularity. From the edge critical exponents of the Fisher edge singularity 
in one dimension we get y\ = 2 and y e h = 2, which also follows from <j e = — | at the Yang- 
Lee edge singularity. The equivalence of the Yang-Lee and Fisher edge singularities is not 
surprising in light of the transformation which maps the temperature and field variables 
into each other in one dimension ||81|| . In two dimensions it is known that cr e = — | from 
which it follows that yf^ — -g-. If we assume y\ = y^ = -jr- for the Fisher edge singularity, 
we have a e = ~, (3 e = — 4, and r y e = I, which are not far from the values of the exponents 
estimated by series expansions. In three dimensions if we take a e = 1.0 then yl = 3, and 
a e ~ (logarithmic singularity), which seems to be the historical trend, then y^ — 3 as 
well. So in d = 1, 2 and 3 we seem to find y| ~ y^. If y\ = y^, the Fisher and Yang-Lee 
edge singularities are not independent behaviors, and they seem to have only one relevant 
variable. 



V. CONCLUSION 



We have investigated the interesting properties of the Q-state Potts model on the simple- 
cubic lattice using the zeros of the exact partition function on a finite lattice and series 
expansions. The critical behavior of the Potts model in the ferromagnetic and antiferromag- 
netic phases has been studied at the same time based on the distribution of the partition 
function zeros in the complex temperature plane. The distribution of the zeros in the com- 
plex temperature plane reveals the possibility of the intermediate phase transitions below 
the order-disorder phase transition of the Potts antiferromagnet. The behavior at the CT 
singularities, where some thermodynamic quantities diverge with the characteristic CT ex- 
ponents acT, Pct, and 7ct, has been discussed using the partition function zeros and the 
low-temperature series expansion of the Potts model in the absence of an external mag- 
netic field. In the complex temperature plane the CT singularities coexist with the physical 
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critical points for H q = 0. 

We have also studied the veiled properties of the Potts model for nonzero magnetic 
field using the partition function zeros and the high-field, low-temperature series expansion. 
As H q increases, all the zeros in the complex Y = e^ J plane move away from the origin. 
For H q > the physical critical points Yf M and Y^ FM disappear and the Fisher edge 
singularities Y e appear in the complex temperature plane. For different values of X = e^ Hq we 
have calculated the locations and the edge critical exponents of the Fisher edge singularities 
of the Ising and three-state Potts models on the simple-cubic lattice. The values of the 
edge critical exponents a e , /3 e , and 7 e at a Fisher edge singularity are universal, that is, 
independent of X and Q but dependent on dimension d. We found that the number of 
pairs of the Fisher edge singularities is equal to d. We have considered the Yang-Lee edge 
singularity X e , which appears in the complex X plane for Y > Y c instead of the physical 
critical point X c = 1, and its edge exponent a e = 1/S e . In one, two, and three dimensions 
the Yang-Lee edge singularity is related to the Fisher edge singularity. They seem to have 
only one relevant variable yf ~ y^. 
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TABLE I: The density of states Qq = y(E) of the seven-state Potts model on the 3x3x3 simple-cubic 
lattice. 

£ tt Q =7(E) 

7 

1 

2 

3 

4 

5 756 

6 378 

61 8297099724190402024200 

62 8931912475351510844640 

63 8562381313020489303720 

70 60500919374783049960 

71 9822183432973269120 

72 794516261799074640 



20 



TABLE II: The ferromagnetic critical temperatures Yf M of the Q-state Potts models on the 
simple-cubic lattice. l^ M (numer.) are the values of the critical temperatures estimated by various 
numerical methods in the literature. 



Q_ 
1 

2 
3 



Y C FM = 1/(Q + V2Q^T) 1/3 



Y 



FM 



V(t _1 + Q 1/3 ) 



Y ( F M (mimer.) 



0.793700. 
0.644689. 
0.575879. 



0.531886. 



l 

2 

0.475240. 



0.455151. 



0.758546. 
0.633620. 
0.568001. 



0.524738. 



0.493027. 
0.468289. 

0.448181. 



0.753[||, || 



0.64190, H| 
0.571 \M\ 



0.5766j57|] 



0.5767@1 
0.5230 
0.524[||] 
0.532 
0.501 
0.472@] 
0.4770 
0.4560 



TABLE III: The locations (Act = Y CT ) and the corresponding critical exponents of the complex- 



temperature (CT) singularities for the three-state Potts models on the simple-cubic lattice. 



Act 


acT 


PCT 


1CT 


a+2/3+7 


0.1036(49) ± 2.0157(35)i 


1.107(8) 


-0.063(1) 


1.161(38) 


2.142(39) 


-1.6158(119) ±0.7958(142)i 


1.208(212) 


-0.063(2) 


1.237(166) 


2.319(269) 
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TABLE IV: The locations and the edge critical exponents of the Fisher edge singularities for the 



Ising (Q 


= 2) 


and three-state Potts (Q 


— 3) models 


on the simple- 


•cubic lattice for 


dull external 


magnetic field X = e@ Hq . 










X 




Ye 


a e 


Pe 


7e 


a e +2/? e +7 e 


100(Q = 




1.249(9) ± 0.642(4)* 


1.079(54) 


-0.418(27) 


1.024(10) 


1.267(67) 






±1.323(12)*' 


1.047(64) 


-0.410(86) 


1.075(78) 


1.303(157) 






-1.249(9) ± 0.642(4)*' 


1.079(54) 


-0.418(27) 


1.024(10) 


1.267(67) 


200(Q = 




1.433(10) ± 0.737(4)*' 


1.046(62) 


-0.422(24) 


1.031(9) 


1.233(71) 






±1.500(12)*' 


1.052(2) 


-0.415(107) 


1.058(125) 


1.280(197) 






-1.433(10) ± 0.737(4)*' 


1.046(62) 


-0.422(24) 


1.031(9) 


1.233(71) 


500(Q = 


2) 


1.649(11) ± 0.876(5)*' 


1.057(24) 


-0.425(24) 


1.038(9) 


1.246(42) 






±1.771(21)*' 


1.050(5) 


-0.418(76) 


1.058(114) 


1.271(157) 






-1.649(11) ± 0.876(5)*' 


1.057(24) 


-0.425(24) 


1.038(9) 


1.246(42) 


100(Q = 




1.217(35) ± 0.573(21)*' 


1.065(141) 


-0.611(50) 


1.003(33) 


0.846(161) 






0.052(13) ± 1.201(32)z 


0.999(155) 


-0.608(14) 


1.076(134) 


0.858(206) 






-1.070(25) ±0.603(17)? 


1.070(158) 


-0.603(17) 


1.064(60) 


0.927(171) 


200(Q = 


3) 


1.342(36) ± 0.660(22)*' 


1.075(35) 


-0.605(16) 


0.998(82) 


0.863(92) 






0.054(11) ± 1.359(38)* 




— U.DlU^ZD J 


l.Ut)4^1o ( J 


n QQOi'oeQ N i 
U.ooy^zoo J 






-1.199(29) ±0.682(21)* 


1.039(227) 


-0.610(20) 


1.066(66) 


0.885(238) 


500(Q = 


3) 


1.535(39) ± 0.787(24)*' 


1.071(17) 


-0.605(48) 


1.012(22) 


0.872(73) 






0.056(8) ± 1.601(41)* 


1.058(134) 


-0.615(32) 


1.076(106) 


0.905(177) 






-1.397(33) ±0.802(22)* 


1.096(66) 


-0.613(37) 


1.063(155) 


0.934(177) 
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TABLE V: The locations and the edge critical exponents of the Fisher edge singularities for the 
Q-state Potts models on the square lattice for X = 100. 



Q 


Ye 


a e 




7e 


a e +2/? e +7 e 


3 


1.232(2) ± 1.048(3)* 


1.217(19) 


-0.197(6) 


1.205(33) 


2.029(39) 




-1.063(13) ± 1.068(6)* 


1.238(26) 


-0.172(25) 


1.257(69) 


2.151(82) 


4 


1.159(5) ± 0.929(10)* 


1.181(76) 


-0.180(4) 


1.121(77) 


1.941(76) 




-0.908(13) ± 0.959(10)* 


1.136(17) 


-0.127(11) 


1.206(16) 


2.088(28) 


5 


1.103(4) ± 0.859(14)* 


1.181(109) 


-0.173(4) 


1.110(98) 


1.946(109) 




-0.834(18) ± 0.889(27)* 


1.131(176) 


-0.109(28) 


1.207(166) 


2.120(245) 
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FIG. 1: Partition function zeros in the complex Y = e~@ J plane of the Q-state Potts models on 
the 3x3x3 simple-cubic lattice. The solid squares (Y c ) show the locations of the ferromagnetic 
critical points estimated in the literature for Q = 2p4|, [55| , 4[6C], and 7 [pi]]. 
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FIG. 2: Partition function zeros in the complex A = e /3J plane of the three-state Potts model on 
the 3 x 3 x 30 simple-cubic lattice. The open triangles show the locations of the ferromagnetic 
(A FM = 1.7342)0 and antiferromagnetic (A^ FM = 0.4424) g(J critical points estimated in the 
literature. The X marks are the suggested locations of additional transition temperatures at 
A KH = 0.3679||, Q, A RL1 = 0.2806 and A RL2 = 0.2736fll]] (Because A RLi and A RL2 are too 
close, they are not clearly distinguished in the Figure. ), and A KRL = 0.22980] below A^ FM 
for the antiferromagnetic model. The open circles are the locations of two pairs of the complex- 
temperature (CT) singularities estimated from the series analysis. 
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FIG. 3: Partition function zeros in the complex X = e@ Hq plane of the three-state Potts model on 
the 3 x 3 x 3 simple-cubic lattice for Y = e~P J = 0.1, Y = Y c « 0.5766 Y = °- 7 > and Y = °- 9 - 
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FIG. 4: Partition function zeros in the complex Y = e~@ J plane of the three-state Potts model on 
the 3x3x3 simple-cubic lattice for an external magnetic field X = e@ Hq . The closed circles show 
the locations of three pairs of the Fisher edge singularities estimated from the series analysis for 
X = 100. 
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